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The full distribution of the conductance P(G) in quasi-one-dimensional wires with rough sur- 
faces is analyzed from the diffusive to the localization regime. In the crossover region, where the 
statistics is dominated by only one or two eigenchannels, the numerically obtained P{G) is found to 
be independent of the details of the system with the average conductance (G) as the only scaling 
parameter. For (G) < e^//i, P{G) is given by an essentially "one-sided" log-normal distribution. 
In contrast, for e^ /h < (G) < 2e^ /h, the shape of P{G) remarkable agrees with those predicted by 
random matrix theory for two fluctuating transmission eigenchannels. 
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The large conductance fluctuations occurring in meso- 
scopic samples have attracted large attention during the 
last decade ||g]. The magnitude of the fluctuations 
{^ e^/h) is almost independent of the mean value of the 
conductance (G) and the system size ('universal conduc- 
tance fluctuations') Q. This holds for wires or waveg- 
uides in the diffusive regime where the length of the 
wire L is much larger than the mean free path i but 
still much smaller than the localization length ^. As the 
length of wire L approaches the localization length (,, the 
conductance fluctuations become of the same order as 
the averaged conductance. The averaged values are then 
not enough to describe the transport properties and the 
knowledge of the conductance probability distribution is 
of primary interest. 

While for a single-mode wire, where £, ^ i, the proba- 
bility distribution P{G) can be computed for any system 
length, it is remarkably difficult to extend this result to 
the N-mode case [0Q|. Based on numerical simulations 
as well as on perturbation theory |]l],||j], PiG) is ex- 
pected to evolve from a Gaussian (deep in the diffusive 
regime) to a lognormal distribution (deep in the local- 
ization regime). However, the behavior of P{G) in the 
crossover between diffusive and localized regimes is not 
well known. Recently, both experimental H and numeri- 
cal p-[l(]|] studies have reported that the conductance dis- 
tribution at an integer quantum hall transition presents 
a striking shape in the form of a broad distribution. Re- 
lated to that, recent numerical studies suggested a uni- 
versal distribution at the metal- insulator transition [0 . 
Based on the Dorokhov-McUo-Pereyra-Kumar (DMPK) 
equation p2| for quasi-lD wires in absence of time rever- 
sal symmetry it has been found [Ol that the crossover 
region between metallic and insulating regimes is highly 
non-trivial, and shows one-sided log-normal distribution 
for the conductance at the transition. It was suggested 
P,13l that these distributions might not be neither a pri- 



vate issue of the quantum Hall transition nor a product 
of the removal of the time reversal symmetry due to the 
presence of magnetic fields. Moreover, numerical calcu- 
lations p in 2D and 3D suggest that the form of P{G) 
at the critical point is independent of the dimensionality 
of the system and of the model. 

Here we present the results of extensive numerical 
calculations of the conductance distribution P{G) and 
transport eigenchannels all the way from the diffusive to 
the localization regime, paying special attention to the 
regime crossover. We show that close to the crossover 
regime, where the averaged conductance (G) is between 
1/2 and 1 (in units of 2e^//i), the general characteristics 
of P(G) are dominated just by two fluctuating eigenchan- 
nels. The diffusion-localization transition is fully charac- 
terized by a transition from two to only one eigenchannel, 
occurring at a critical value {G)c ~ 1/2. Right at the crit- 
ical value (G)c, the distribution is almost perfectly flat 
and almost identical to that obtained for 2D systems in 
the quantum Hall regime [B|. In the insulating regime 
((G) < {G)c) our results are consistent with the one-side 
log-normal distributions predicted by the DMPK equa- 
tion |l|]. 

Most of the previous work on transport properties of 
disordered wires has been focused on the influence of mul- 
tiple scattering with bulk defects. In this work we are 
going to deal with surface corrugated wires (SCW) as a 
model system. The interest on these systems is twofold. 
On one hand, for very thin wires, the main source of mul- 
tiple scattering comes only from the surface roughness. 
On the other hand, in absence of bulk defects, the trans- 
port through SCW, like in the DMPK approach, does 
not contain the effects of wave function correlations in 
the transverse direction. Many statistical properties of 
transport through SCW |13-M] have been found to be 
in good agreement with those predicted by random ma- 
trix theory (RMT) for bulk defects. In particular, the 



analytical distributions of transmission coefRcients | [l9| ] 
obtained from DMPK equation are nicely reproduced in 
numerical simulations of SCW |l6[ . Our model system 
consists of a surface corrugated two-dimensional (2D) 
wire (see inset in Fig. 0). The corrugated part of the 
wire, of total length L and perfectly reflecting walls, is 
composed of slices of length I. The width of each slice has 
random values uniformly distributed between Wq — 5 and 
Wq + 5 about a mean value Wq. The transmission ma- 
trix t is exactly calculated for each realization by solving 
the 2D wave equation by mode matching at each slice, 
together with a generalized scattering matrix technique 
|^,Q. The conductance is given by G = tracejtt^}. To 
obtain the mean values ((. . .}) and the probability dis- 
tributions we have used 10000 independent realizations 
of the disordered guide. 

We have performed calculations for three different val- 
ues of Wo/A: 4.9, 2.6 and 1.8, allowing TV = 9, 5 and 3 
propagating modes respectively (A would be the Fermi 
wavelength for electron transport or the wavelength of 
the diffusive incoming waves in the case of transport of 
classical waves). We have considered 1/5 = 3/2 and dif- 
ferent ratios Wq/5 = 7 (for TV = 5, 3) and Wq/S = 13.25 
(for A^ = 9). The behavior of different mean values of 
transport coefficients as a function of the length L has 
been discussed previously in detail |l^Jl^. In the dif- 
fusive regime, the averaged resistance follows a typical 
ohmic behavior ((1/G) ~ 1/N + L/£^) while, in the local- 
ization regime, (InG) ~ —L/^. In Figure 1 we plot our 
results both for (G) (Fig. la) and for (In G) (Fig. lb) 
versus s = NL/S^ for TV = 9 (dashed line), N = 5 (solid 
line) and TV = 3 (long-dashed line). The corresponding 
localization lengths are C = 187 Wo (iV = 9), ^ = 34.4Wo 
(TV = 5) and ^ = 51Wo [N = 3). 

The analysis of the conductance distributions P{G) 
along the transition from the diffusive regime L < ^ to 
the localization regime L > S, (see Fig. ||) shows that 
the distributions are far from evolve smoothly from the 
Gaussian to the lognormal distribution as one could ex- 
pect from the behavior of the transmission coefficients 
|^,|6[. Our results for P{G) are shown in Figure || for 
some selected values of (G) around the transition from 
the diffusive to the localization regime ((G) w 1 (a), 4/5 
(b), 1/2 (c) and 1/3 (d); those values correspond to the 
symbols in Fig. 1). The results for iV = 9 (O), iV = 5 (o) 
and A'^ = 3 (■) are almost identical (within our numerical 
accuracy) suggesting that the only parameter determin- 
ing the shape of P{G) is the average conductance (G) as 
expected from scaling arguments p3]. 

It is remarkable that the almost perfectly flat distri- 
bution obtained for (G) ~ 1/2 (Fig. 2c) is in full agree- 
ment with those obtained for 2D systems in the quan- 
tum Hall regime at the critical point P,M: a flat dis- 
tribution except for a small dip at G « [pilj] and 
an exponential cutoff for G > 1. As the averaged con- 
ductance goes slightly below 1/2, P{G) becomes a "one- 



sided" log-normal distribution with a sharp cut off at 
(G) > 1 (Fig. 3) in agreement with both numerical [0J|] 
and analytical [|3[ results. In Figure 3 we have plot- 
ted our results for P(lnG) for different values of (G) to- 
gether with the best fit to a one-sided log-normal dis- 
tribution (thick solid lines). Although the agreement 
is very good for all values of (G) < 1/2, the distri- 
bution at the critical value of 1/2 is clearly better de- 
scribed by P(ln G,(G) = 1/2) = exp(lnG) (thin line 
in Fig. 3a) which corresponds to a fiat distribution 
(P(G, (G) = 1/2) = 1 , < G < 1). Although there 
is no phase transition in quasi-one-dimensional systems 
(ID), our results suggest a special generic behavior at a 
critical {G)c = 1/2 which is certainly not consistent with 
a one-sided log-normal distribution. 

In order to get more insight on the exhibited behavior 
of the conductance distributions we have made an exten- 
sive analysis in terms of the transmission eigenchannels 
JTi}. These eigenchannels are the eigenvectors of the ma- 
trix tt^ and they form the natural basis to analyze the 
properties of the conductance (G = trace(tt^) = J^i '''«)• 
Following Imry's work pl| , we arrange the eigenvalues is 
descending order, and perform the averaging only taking 
those eigenvalues of the same level into account (e.g. (ti) 
corresponds to the average of the highest eigenvalue of 
each realization). In Fig. 4 we show the evolution of the 
average value of the eigenchannels as a function of (G) . 
In the diffusive regime ((G) > 1), in agreement with the 
behavior discussed by Imry, almost all the eigenchannels 
are either fully open or closed and only a few of them fluc- 
tuate giving rise to UCF. However, their exact behavior 
depends on the particular parameters of the wire. Close 
to the crossover to localization, when all the eigenchan- 
nels are closed and only two (or one in the localization 
threshold) of them are partially open, the average eigen- 
channel transmissions (Fig. 4) as well as the conductance 
distributions (Fig. 2) do not depend on the size of the 
wire or the defect details. In this sense, our results sug- 
gest an universal behavior of the conductance with (G) 
as the only scaling parameter. 

In contrast with the insulating regime (G) < 1/2, 
where the conductance is known to be well described by 
one-sided log-normal distributions U^, in the crossover 
regime (G) > 1/2, where the statistics is dominated by 
one or two fluctuating channels, there is no analytical 
result available for the conductance distribution. As a 
first analytical approach to this problem, given our lim- 
ited knowledge of the statistical correlations between dif- 
ferent eigenchannels, a possible choice of the statistical 
ensemble is that which maximizes the information en- 
tropy subject to the known constrains of flux conserva- 
tion and time-reversal invariance. In the RMT context, 
this leads to the circular orthogonal ensemble. In our 
case, we also know that, close to the onset of localization 
(for 1/2 < (G) < 1), only two eigenchannels actually con- 
tribute to the conductance, independently of the initial 



number of channels N. From the polar decomposition of 
the scattering matrix [p|p4| it is possible to write the joint 
probability distribution of n transmission eigenchannels 
(with Ti — 0, i ~ n+ 1, . . . , N) as 
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where a = (1 + n)((G) - n/2)/{n - (G)) ^. Within 
this approach, the conductance distribution is fully de- 
scribed by the number of fluctuating channels n and the 
mean value of the conductance (G) [^. Simple closed 
expressions for P{G; {G),n) can be obtained for n = 1,2, 
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In Fig. g we have plotted the RMT distributions (con- 
tinuous thick line) for different (G) {n — 1 for (G) < 1/2 
and n = 2 for 1/2 < (G)), together with our numerical 
results. As it can be seen, there is a very good agree- 
ment between the analytical results and the numerical 
calculations. In particular, our RMT approach captures 
some important features of the distributions such as the 
almost perfectly flat distribution at (G) = 1/2 and the 
cusp point in P{G, (G)) at G = 1. Since the RMT re- 
sults are quite general, we expect qualitative similar con- 
ductance distributions in higher dimensions close to the 
critical regime. 

In summary, we have analyzed the evolution of the 
conductance distribution all the way from the diffusive 
to the localization regime. Close to the crossover regime 
we have shown that the distributions are independent of 
the details of the system with the averaged conductance 
(G) as the only scaling parameter. For (G) between 1/2 
and 1, the statistics is dominated by one or two fluctuat- 
ing eigenchannels and the numerical conductance distri- 
butions are surprisingly well described by RMT results 
At a critical value {G)c = 1/2, the distribution is almost 
perfectly flat. In the insulating regime ((G) < 1/2) our 
results are consistent with one-sided log-normal distri- 
butions. The similarities between our results and those 
obtained in very different situations suggest that the con- 
ductance distributions exhibit a universal behavior at the 
crossover regime. 
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FIG. 1. Averaged conductance (G) (a) and averaged log- 
arithm of the conductance (InG) (b) as a function of the 
length of the disordered part of the wire for VKo/A = 4.9, i.e. 
9 modes in the clean part (dashed line), Wo/X — 2.6, i.e. 5 
modes (solid line) and Wo/A = 1.8, i.e. 3 modes (long-dashed 
line). Vertical lines indicate the localization length for each 
case. Symbols are the points at which the conductance dis- 
tributions have been calculated. Inset: Schematic view of the 
system under consideration. 




FIG. 2. Conductance distribution for the different averaged 
conductance values [(a){G) = 1; (b){G) = 4/5;(c){G} = 1/2; 
(d){G) = 1/3] shown in Fig. Q. Symbols are for iV = 9 (O), 
N = 5 (o) and A'' = 3 (■). Continuous lines represent the 
analytical results of our random matrix approach. 
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FIG. 3. Distribution of InG for different values of the 
average conductance beyond the critical point {{Gc} ^ 1/2). 
Thick solid lines are the best fits to a one-sided log-normal 
distribution. Thin solid line in (a) corresponds to 

P(lnG, (G) = 1/2) = exp(lnG) (i.e. P{G,{G) = 1/2) is a 
uniform distribution). 
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FIG. 4. Averaged transmission eigenvalues (ji) as a func- 
tion of the averaged conductance. Dashed line corresponds to 
a 9-mode wire, solid line to a 5-mode wire and long-dashed 
line to a 3-mode wire. 



